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Abstract 

We present an analytical study of the Lyapunov exponent A for electronic 
states in a disordered tight-binding ring, with sites of spacing c, in the 
presence of a non-hermitian field component h, for weak disorder. This system 
has been proposed by Hatano and Nelson (HN) as a model for the depinning 
of flux lines bound to columnar defects in a superconducting cylindrical shell 
in a transverse magnetic field proportional to h. We calculate the Lyapunov 
exponent to second order in the disorder in the domain of complex energies 
in an ordered ring in a non-hermitian field. The result exhibits an important 
symmetry property, A(— /i) = — A(/i), which reflects the directionality of the 
non-hermitian localization. By combining the results for A with a previous 
calculation of the complex eigenvalues for weak disorder, we show that the 
localization length of the complex eigenenergy states is generally larger than 
the circumference (Nc) of the ring. This supports the initial suggestion of 
HN that these eigenstates are delocalized. Our perturbative treatment of 
the disorder when applied for an open chain embedded in an infinite non- 
disordered chain recovers the Thouless formula for the zero-field localization 
length. 
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I. INTRODUCTION 

It has been known since the classic papers of Mott and Twose||l]| and of 
BorlandQ that disorder localizes electronic states in an open linear chain. This 
observation was reinforced about two decades later, by detailed analytical calcu- 
lations of the inverse localization length (ILL) for weakly disordered tight-binding 
open chains [3-7], for energies in the band of the ordered system. 

Recently, renewed interest in one-dimensional localization has arisen in two 
important contexts involving closed circular chains (rings) rather than open linear 
chains. The first one is the problem of the universal persistent current induced in 
thin mesoscopic metallic rings threaded by an Aharonov-Bohm flux, which has been 
studied both theoretically [8-10] and experimentally [|l^. The other one relates 
to flux lines pinned to random columnar defects in a superconducting cylindrical 



shell[[12[. This system can be modelled by a non hermitian quantum hamiltonian 
for an electron in a disordered ring in the presence of a constant imaginary vector 
potential (IVP) (proportional to the external magnetic held acting transversely to 
the flux lines in the superconductor||l^). The localization-delocalization transition 
(corresponding to a pinning-depinning transition for the flux lines) predicted by 



Hatano and Nelson |12|] at a threshold vector potential, has generated considerable 



interest from an early stage[[T^, [Tj[ till recently [[T5|, |T6|. Localization in the presence 
of an imaginary vector potential is referred to speciflcally as directed localization[0 
because, unlike in a real vector potential, the symmetry between right and left 
moving particles is broken while time reversal symmetry is not. 

The previous analyses of the localization-delocalization transition in the non- 
hermitian hamiltonian of Hatano and Nelson[12-15] are mainly based on the study 
of the eigenvalues which may be complex. Indeed it has been argued that real 
eigenvalues correspond to localized states while complex eigenvalues correspond to 
extended states[0]. In particular, both the analyses of the directed-localization 
problem and those of the effect of disorder on persistent currents [9,10] use the fact 
that in zero field the electronic states in a one-dimensional disordered system are 
localized and correspond to real eigenvalues. 

The purpose of this paper is to discuss the theory of localization and of 
directed localization by means of a first principles analytical study of the Lyapunov 
exponent A of the electronic states in a weakly disordered tight-binding ring in a 
non-hermitian field. The ILL in the non-hermitian field, which we shall refer to 
more specifically as the inverse directed localization length (IDLL), 1/^, is related 
to A by 1/^ = |A|. 

The Schrsdinger equation for the hamiltonian of a disordered one-orbital 
tight-binding ring of circumference Nc (with N the number of sites and c the lattice 
spacing) in an IVP reduces to the set of equations 

e'^Vn+l + e-'^Vn-l + Sn^n = E(pn , U = 2, 3, . . . N - 1, (1) 

e'"'ip2 + e-'^Viv + eiifi = Eifi, (2) 

e'^Vi + e"''Vjv-i + Sn^n = E^^- (3) 

Here ipn is the amplitude of a wavefunction at site n, E and the parameters Sn are 
the energy and the random site energies in units of minus the nearest-neighbour 
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hopping parameter, respectively. The hamiltonian in is non hermitian because 
of the real field parameter h, corresponding to an IVP||12||. 

In an interesting recent paper, Brouwer et al. presented an analytic study of 
the complex eigenvalues of the non-hermitian tight-binding ring for weak disorder 
[14]. This study is complemented here by a corresponding analytic study of the IDLL 
for complex energy states for weak disorder. Our treatment is valid for sufficiently 
strong fields such that the disorder may be treated perturbatively. The localization 
length of complex energy states as defined in a ring of infinite circumference 
[Nc = oo) is generally finite. This is similar to the decay length ^ of the transmission 
coefficient introduced by Brouwer et al. in their analysis of complex eigenvalues 
which is also found to be finite and given by ^ = (2|/i|)~^[0]. However, when 
applying these results to a finite ring one must consider the restriction on the validity 
of the perturbation treatment of the disorder which is imposed by the discreteness 
of the energy spectrum of Specifically, the corrections to the zeroth order 

complex energy levels due to the disorder must be small compared to the level 
spacing of the ring. Using a previous calculation of the eigenvalues of (|l]-^ for 



weak disorder, this condition allows us to show that the localization lengths of the 
complex eigenenergy states are generally larger than the circumference of the ring 
at the considered field strengths. These results explicitly support the suggestion of 
HN [0 that the complex energy states in a finite ring are delocalized. 



The paper is organized as follows. In Sect. II we present the detailed solutions 
for the amplitudes ipn in for weak disorder, at arbitrary energy E both in the 

domain of the complex spectrum of the ordered system, and in a range of real 
energies. From these solutions we obtain the explicit expressions for the Lyapunov 
exponent for the eigenstates of the ring, for /i 7^ 0, to second order in the random site 
energies. In Sect. Ill we first apply a similar perturbation analysis for weak disorder 
to rederive the Thouless formula for the zero field ILL for the case of an open chain. 
This further supports the validity of our calculation for the case of a ring. Next we 
discuss the properties of our general expression for the IDLL and its application to 
complex eigenenergy states in finite rings. In Sect. IV we conclude with some final 
remarks. 

11. THE NON-HERMITIAN LYAPUNOV 
EXPONENT FOR WEAK DISORDER 

In solving (|l|-0) recursively, initiating the recursion at site n = 1 with an 
arbitrary (fi, it is convenient to define the amplitude ratios i?„ = (pn/^n-i , n = 
2,3, ... N, and Qi = fi/ifN obeying the equations 

e'^^Rn+i + e-'^'R-^ = E~en, ,n = 2,3, . . . N - 1, (4) 
e'"'R2 + e-'"'Q:,^ = E-ei, (5) 

ef^-Q, + e-^'R],' =E-eM, (6) 

which determine the quantities Rn , n = 2,3, ... N, and Qi as function of an 
arbitrary energy E, if one disregards the secular equation giving the eigenvalues, 
which may be written in the form 



2 



N 

e'^'R^iE) + e-'^^ J] R^{E) = E - e,, (7) 

m=2 

using (|2|) and the general definition 



<^n = n ^-'^i , n = 2,3,...iV . (8) 



m=2 



The secular equation (0) will be automatically verified by the quantities Rn{E) 
above whenever E coincides with one of the eigenvalues E = E^. 

The reason for studying the exponential variation of the quantities and 
of the corresponding amplitudes as a function of a generic energy E^ at sites 
n further and further away from the initial site, is that it permits to obtain the 
Lyapunov exponent of the electronic eigenstates of the ring on the basis of the Mott- 
Twose-Borland conjecture||I], |18|. Indeed, Mott, Twose and Borland argue that 



when E is close to an eigenvalue the exponential rate of variation of the amplitudes 
(Pn at far away sites n tends to the (largest) exponential rate of localization of the 
corresponding eigenstate of the ring about a fixed localization centre. The Lyapunov 
exponent which characterizes the exponential growth (or decay) of at distant sites 
(nc — » oo), at real or complex energies E, is given by |[T^] 

A = lim — In 

{n<A')^oo nc 



1 " 

= lim — Vln|i?„| . (9) 

n^oo nc 

and since A self-averages to a central limit value||2n| we further have 



1 " 

A = lim — y(ln|i?p|) . (10) 

n-*oo nc 

p=2 

The Lyapunov exponent may generally take positive- as well as negative values. 
Since all sites of a ring are equivalent by definition, there is no fundamental dis- 
tinction between states corresponding to different signs of the Lyapunov exponents. 
In contrast, for an open chain, a positive Lyapunov exponent is associated with 
states which are growing as one is moving away from an edge site n = 1 while a 
negative Lyapunov exponent corresponds to a state which decays with increasing 
distance from the edge. The former state is generally referred to as a localized state 
while the latter might be termed an antilocalized state||2l| to emphasize the fact its 
amplitude is smallest deep in the interior of the chain. 

Assuming the disorder to be weak, we now analyze the perturbative solutions 
of (4-6) to second order in the site energies Si. Since we are interested in the study of 
the localization length, we only consider solutions as a function of a generic complex 
energy E, thus ignoring the consistency-or eigenvalue equation (|^). In the absence 
of disorder the solutions of the non-hermitian system (4-6) are of the form 
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=g; = e^^ n = 2,3...iV, (11) 
and correspond to a generic complex energy 

E = 2cos{q-ihc) (12) 
expressed in terms of an arbitrary real wavenumber q in units of the lattice para- 



meter. The equation ([121) defines energies lying on an ellipse, 



{ReEf ^ {ImEf _ ^^^^ 
cosh^ he sinh^ he ' 

whose semi-axes 2cosh/;,c and 2sinh|/i|c give the half- width of the energy band 
of edges ±2 cosh he on the real axis and of the corresponding band with edges 
±sinh|/i|c on the imaginary energy axis, respectively. The meaning of the ener- 
gies ([l2|) is that they correspond to the domain of extended Bloch state eigenvalues 
in the absence of disorder, which are given by the special values for q = 2'7ik/N , 
/c = 0,1,2,... as a consequence of single- valuedness of the wave functions on the 
ring, expressed by (|^. On the other hand, the meaning of ImE has been discussed 
by Hatano and Nelson in the context of depinned flux lines [p!^ . In fact, Hatano 
and Nelson argue convincingly that complex energies are a necessary feature of the 
delocalized states which may appear in the depinning transition. 

The corrections in i?„ and Qi at first and the second orders in the energies 
ei,i = 1^2,...N, are determined by inserting the expansions 
Rn = e"^ + R^n^ + n = 2, 3, . . . AT , Qi = 9 + Q^^ + Qf\ involving lin- 

ear and quadratic terms, in the equations (4-6) and identifying terms of the same 
order in the order by order expansions of them. This leads to the following sets of 
equations at first and second orders, respectively: 



R^^l^ - ~a R^^^ = -e-'^^En , n = 2, 3, . . . - 1, (14.a) 

- i?!v^ = -e-'^'id SN + £i), (14.b) 

- d = -a e-^^i?i^)2 , n = 2, 3, . . . - 1, (15. a) 

- d^ i?Sv^ = -d^{l + d) e-*'?i?S^^ + dVd eN{2d R^^^ - c'^'en). (15.b) 

Note that the quantities Q^^^ and Qf^ have been eliminated since they do not enter 
in the determination of the wavefunction amplitudes (||). In the above equations we 
have defined 

d = exp(— 2ig) , q = q — ihe, (16) 

where g in a wavenumber with an imaginary part arising from the non-hermitian 
field. 
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By solving successively the coupled first and second order diff'erence equations 
1^14. a, b) and (15. a, b) for the quantities Rn^ and Rn^ we get: 



p(i) 



n-1 



m=2 



-he / ^-1 



m=2 



(17.a) 
(17.b) 



m=2 ^ ^ 



R 



(2) 



2 Af 



n = 3,4, ...A^, 



m=2 



(18.a) 
[IS.b) 



It should be noted that the expression for R2'' and R2' , equations (17. b) and 
(18. b), are invalid at zero field, where they diverge at the energies E = 2cosg, 
q = , A; = 0, 1,2, . . . corresponding to the eigenvalues of the ordered ring for 
h = 0. In fact, as is well known, perturbation theory for eigenstate energies in a 
weakly disordered ring also diverges for h = 0^^. This strongly suggests that the 
Lyapunov exponents of the zero field eigenstates in the ring may not be discussed 
by perturbation theory. 



?(2) 



The next step is then to expand ln(i?e Rp) in (|10D up to second order 



terms 22 



\n\Rp\ = ^{e-'ml^^ + c.c] 



1 

+ 2 



,-iq 

_^(1)2 

2 P 



+ c.c. 



p = 2,3,...N, 



(19) 



and to obtain the averages of the quantities on the r.h.s. using the explicit ex- 
pressions (17. a), (17. b) and (18. a), (18. b). Here we assume the 5j0s to be identi- 
cally distributed, independent gaussian variables with mean zero and a correlation 
(eiEj) = £q 6ij, where the choice el = would correspond to using the familiar 

rectangular distribution of width W. The determination of the non-zero averages 

(1)2 (1)2 

proceeds by first finding the averages {R\' ) and {Rj;'^) , p^2, from (17. a) and 

(2) 

(17. b), which are then inserted in (18. a) and (18. b) for obtaining the averages {R2 ) 

(2^ 

and ) , p ^ 2. The sums over lattice sites in these various averages are geomet- 
ric progressions which are easily handled. The final form of the averages of interest 
are remarkably simple since we obtain 



(2)\ 



4 e 



-2hc 



1 + a 



~N 



;i-a)(l-a2) Vl-a^ 



, P = 3,4, 



(20.a) 

N-l. 
(20.b) 
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This reflects the independence of these quantities of the initial site(s) for the recursive 
solution of the tight-binding equations, as expected. By substituting these results 
in (|T9| ) and ( |T0| ) we obtain our final expression for Lyapunov exponent (A^ oo) to 



second order in the disorder, 



4 (1 - a)2 \l~a 



+ C.C., (21) 



where we have used the definition ([Tq) . The equation (|2lD is valid for complex 
energies corresponding to the energy band ([T3|) of the pure ring in a non-hermitian 
field h^O. 

The above analysis may be extended in principle to the description of the 
IDLL for real eigenvalue states. It suffices to choose the parameter q in (|TlD to be 
complex with the forms 



q = in , (22. a) 

q = 7T + iK , (22. b) 

where k is real. Indeed the values (22. a, b) define generic real energies 

E = ±2 cosh (k - he) , (23) 



which together with ([TTD and (22. a, b) are solutions of (4-6) for the pure system. 
Eigenvalues of the non-disordered ring obtained by inserting ([TT|), (22. a, b) and ( pSj ) 
in (1^ correspond to k = i.e. to the real band edges ±2 cosh he. On the other 
hand, the perturbation of eigenstate energies for weak disorder does not lead to real 
eigenvalues apart from the edges (±2 cosh he) of the Re i?-domain of the complex 
eigenvalues [14,17]. Therefore the perturbation study of the IDLL for real energies 
is of no interest. 



III. DISCUSSION OF RESULTS 

A. Validity of perturbation theory 

In this subsection we demonstrate the validity of the perturbation expan- 
sion of Sect. II by showing that when applied to the case of an open disordered chain 
it recovers the well-known Thouless formula for the inverse localization length in 
the absence of an IVP. Here an open disordered chain is defined as a segment of 
length Ne which is connected at both ends to semi-infinite non-disordered discrete 
chains. The presence of semi-infinite non-disordered chains is indeed required for 
the existence of a transmission coefficient T related to the inverse localization length 
by T f« e~^^'^ for oo||2^. This implies that for an open chain the solution of 



equation (|J) in the absence of disorder is still given by (for h = 0) 

= E = 2 cosg , n = 1, 2, . . . AT. (24) 
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In order to study the Lyapunov exponent for a large open chain for weak disorder 
it is then sufficient to solve the perturbation equations (14. a) and (15. a) alone in 
terms of arbitrary finite ipi and R^2 \ ^2'^ fhe expansion 

i?2 = + R^^'> + i?f + . . . . (25) 



The relevant averages in ([lOD , using ([191) are obtained from (17. a) and (18. a) for 
/i = and are given by (with a = e~^*^) 



(26) 



+ — ^_^^ |l+a--a>-^-.-^n (27) 



By inserting these results in ([19|) and noting that only terms independent of p 
contribute to the Lyapunov ( [T0| ) we finally obtain, using (|24 



^4 24^ • (^«) 

which is the well-known Thouless formula^, for the ILL showing that all eigen- 
states with eigenvalues lying in the band (^) are exponentially localized (with the 
usual proviso that the Mott-Twose-Borland argument about the connection between 
eigenstates and exponentially growing solutions from both ends of a long open chain, 
at an energy E, is valid). 

B. Directed localization on a ring 

In the presence of an IVP, the rates of hopping from a given site to its nearest 
neighbour to the left- and to the right hand sides on the ring are different. This 
asymmetry of the hopping rates is reflected in the antisymmetry of the Lyapunov 
exponent of the eigenstates under reversal of the non-hermitian field h, 



X(-h) = -X{h), (29) 



which is demonstrated explicitly by the equation (^), to lowest non- vanishing order 
in the gaussian disorder. To our knowledge this important antisymmetry property 
has not been discussed before. The physical interpretation of ([29|) is simple. We 
might expect the Lyapunov exponent in ( pl| ) to describe a symmetric exponentially 
localized state: this means that since \{h) determines the exponential growth (if 
X{h) > 0) of an initial amplitude ipi with increasing n = 1,2,3, .. . then the cor- 
responding exponent describing the growth of the localized state with decreasing 
n = N, N — 1, N — 2, . . . should be —X{h). The fact that the Lyapunov exponent 
defining the exponential evolution of the initial wavefunction with decreasing 
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n = N, N — 1, N — 2, . . . is indeed — A(/i) is explicitly verified in Appendix A. On 
the other hand, the Schrsdinger equation shows that the rates of hopping from 
a given site to its left- and right neighbours, respectively, are interchanged when h 
is replaced by —h. This implies that the rates of exponential evolution of the initial 
wavefunction are also interchanged when h is replaced by —h i.e. the antisymmetry 
property (p9|) holds. Finally, we recall that according to the Mott-Borland [1,2,18] 
matching condition the symmetrically growing (or decaying) initial wavefunction in 
the two directions on the ring, respectively, lead to proper localized eigenstates at 
energies corresponding to complex eigenvalues of (^^. 

The further detailed properties of the IDLL defined for an infinitely large ring 
(Nc oo) are as follows. For N\h\c >> 1 ( |2TD reduces to the central limit result 

A = Aoo sign h, (30) 
where Aoo, which is independent of the sign of h, is given by 



el a 



Aoo — — -rjz TT^ + c.c. 

4 1 — ar 



f 2 / 1 



+ C.C. , (31) 



4 V4-^2 

and defines the IDLL 1/^ = |A| = |Aoo| at complex energies E of the form ([121). 

In the following discussion we regard the Lyapunov exponent A and the IDLL 
as functions of ReE alone, which are obtained by eliminating ImE in favour of ReE 
using ([T3[ ). The main properties of the Lyapunov exponent ( ^0| - [3l[) are: 
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i. For h > 0, X is positive in the interval 



-Ea < ReE < Ea, (32) 

where 

/""^""f . (33) 
V 1 + tanh he 

It possesses a maximum of magnitude 



^ = § (34) 

8 1 + sinh he 

at the band centre, ReE = 0, and vanishes at the values ReE = ±Ed. The equa- 
tion (p4|) shows that the IDLL decreases with increasing \h\e within the domain 
N\h\e ^ 1. Finally A is negative in the domains extending from Ed to the upper 
band edge, 2 cosh he, and from —E^ to the lower band edge, —2 cosh he, respectively; 
at the band edges its value is 



(35) 

8 sinh he 



ii. For /i < 0, A is negative in the interval (^), with a minimum of magnitude (^), 
and vanishes at ReE = ±Eci. A is positive in the complementary intervalls of ( ^2]) 
contained in the real energy band —2 cosh he < ReE < 2 cosh he and takes minus 
the value ( pS]) at the band edges. 

It follows from the above analytical results that the eigenstates of an infinite 
ring in the complex energy domain (pTSD are localized, except possible eigenstates 
with eigenvalues ReE = ±Ed, whose localization lengths would be diverging. On 
the other hand, in the limit of very strong fields {\h\e >> 1) where the effect 



of the disorder becomes negligible, the IDLL ([31|) vanishes asymptotically, at the 
considered complex energies, as 



A ^ - ^ cos 2g sign h e-^''' , (36) 

so that the wavefunction amplitudes, = \ipn\ e*^" = e*''" are delocalized in this 
limit [12]. 

Application to finite rings 

We now turn to the application of the above results in the case of a finite ring 
of circumference Ne. Here only states whose localization lengths are shorter than 
Ne are said to be localized while eigenstates with localization lengths (finite or not) 
larger than Ne are referred to as delocalized. The other important aspect is that 
the discreteness of the energy spectrum of a finite ring imposes a specific restriction 
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on the validity of the perturbative study of a weak disorder. This restriction exists 
over and above the general strong field requirement which is necessary in order that 
the disorder may be treated as a small perturbation in both the real and imaginary 
parts of the energy. We use the results of our recent calculation of the eigenvalues of 
([l|-|^) to explicitate the supplementary condition which must be obeyed by the 
perturbation parameter of the disorder in a finite ring. This condition will allow us 
to derive an upper bound for the IDLL of the form (pi]). 

The discrete zeroth order eigenenergies of the ring are given by 



27rk 

-^o(Q'fc) = 2cosh/iccosgfc + 2isinh/icsingfc , gfc = , A; = 0, ±1, . . . . (37) 

For weak random site energies with zero mean value and a gaussian correlation, 
{siEj) = £q Sij, the averaged second order correction to the energy levels (pTl ) is 
given by 



(fl'-)te-)):.4 '=°^^'^^-;'"' . (38) 

4 sm[qk — ihc) 

The standard heuristic analysis of convergence of non-degenerate Schrsdinger per- 
turbation theory |Q now shows that {E^'^\qk)) must be small compared with the 
separation AEo^qk) = -Eo('?fc+i) — Eo{qk) of the unperturbed energy levels (p7|) for 



convergence to be ensured. By separating the real and imaginary parts in (|37| - |38D 
the convergence conditions of the perturbation theory read 



\Re{E^^\qk))\ « \AReEo{qk)\ , (39.a) 
\Im{E^^\qk))\ « \AImEo{q,)\ , (39.b) 



where 



All 

AEo(qk) — ——(— cosh he sin qk + i sinh he cos qk) , 

Ne 

,^(n\, ^ci (— cos qk sinh he + i sin qi: cosh he) 

{E^ >{qk)) = —sign h — 2 

4 sm qk + smh he 



(40) 
(41) 



for N\h\e» 1. 

It follows from (|40|-41) that the inequalities (^.a,b) are not both fulfilled for 
qk = 0,71 and for = 7r/2 i.e. at the band edges along the real and complex energy 
axes, respectively. For illustration of the condition (|^.a,b) at typical intermediate 
gfc-values we choose qk = 7r/4. For this value the Lyapunov exponent given by ( pT]) 



and (12) takes the form 



= SttttT ' = T • (42) 



4cosh^(2/ic) 



Using (41-^21), the inequality (|39| .a) may be rewritten in terms of the localization 
length ^ = |Aoo|~^ in the form 
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A^c I tanh /icI cosh^(2/;,c) vr , , 

^ » — — - , qk = - , 43 

^ 27r 1 + 2 sinh^ he A ' ^ ' 



while p9|.b) takes the same form with |tanh/;,c| replaced by |coth/ic|, thus imply- 
ing a larger lower bound for ^ than (|43D . As mentioned earlier the perturbative 
treatment of the disorder is generally valid for sufficiently large fields only. In this 



context, (|43D shows that the lower bound of ^, which increases with increasing h, is 
already larger than the ring circumference for fields of relatively moderate strength, 
of the order of he = 1.3. This leads to the important conclusion that the complex 
eigenenergy states in a finite ring are delocalized (and correspond to depinning of the 
associated flux line system ||12]) in the domain of field strengths where a perturbative 
treatment of the disorder is valid. These results provide a detailed analytical confir- 
mation of the suggestion of Hatano and Nelson that the complex energy eigenstates 
are delocalized. 



IV. CONCLUDING REMARKS 

We have calculated analytically the Lyapunov exponent, X{h), for the eigen- 
states of a random circular chain in a strong non-hermitian field h, for complex 
energies, to order in the disorder. 

In a finite ring the discrete nature of the energy spectrum entails a specific 
restriction on the validity of the perturbative treatment of the disorder which implies 
the existence of a lower bound for the localization length ^ = |A(/i)|^^ in (^). This 
shows that the localization lengths at intermediate energies within the complex 
energy band are larger than the circumference of the ring at the strong fields of 
interest in the perturbation theory. It follows that the complex energy eigenstates 
in a finite ring are delocalized, in agreement with the suggestion of Hatano and 
Nelson 

The analytic argument of HN for showing the existence of a transition from 
localized to delocalized states in an IVP is based on an imaginary gauge transfor- 
mation, (p{x) = exp{gx/h)'ilj{x) (with g/h = h), which eliminates the IVP, —ig, 
from the Schrsdinger equation for their continuous random hamiltonian for a ring 
of length Lx- In their analysis HN implicitly replace the twisted boundary condi- 
tion (which is well-known in the context of real vector potentials [25,26]) for the 
transformed wavefunction, namely 



ipix + L^) = e 'i^tp{x) , (44) 

by the periodic boundary confition, i/j{x + L^) = ipix), appropriate for a ring in the 
absence of g. In general this approximation is valid only for a very weak IVP such 
that ^ << l/L^. HN have obtained localized right and left eigenstates, for weak 
IVP, from the localized states (with localization lengths ^„ < Ne) of the ring for 
g = 0. These states are insensitive to the actual boundary condition (0) if e.g. the 
localization length in the case of the right eigenstates, (1/^n — \9\/^)~^ HHIj is small 
compared to Ne. In this case therefore the validity of the replacement of (^4|) by 
a periodic boundary condition extends to a range of somewhat larger values of g, 
namely \g\/h « 1/^n > ^/Lx- It remains, however, that the analytic discussion of 



11 



HN is restricted to a domain of weak fields for the real eigenvalue localized states. 
In contrast the study of energy eigenstates in Sect. II and III applies in a domain of 
large fields, \h\ > 1/c, where the detailed nature of delocalized complex eigenenergy 
states of a finite ring has not been studied analytically before. 
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Appendix 



Here we calculate the Lyapunov exponent describing the evolution of the 
wavefunction amplitude at the successive neighbours A^, A^ — 1, A^ — 2,..., A^ — 
n, . . . of the site 1 along the direction opposite to that considered in Sect. II. From 
successive recursions of (14. a) we obtain 



d(1) _ ~-n 



^Si^+e-'^E^"""''^^-™ ' ^=1, 2,... , Ar-2 , (A.l) 



m=l 



and from (14. b) and ([A.l|) taken at n = A^ — 1, 



(1) 

■N 



-he 



Similarly, from (15. a) we find 



7V-2 

a^'^{a£N + £i) + 



~m-N-l^ 



m=l 



(A.2) 



R^N-n = ~^""Rn^ + e-'^ ^""'"^N-m , U = h 2 AT - 2 



-5 -"5 • • • 5 



m=l 



and from (15.b) with R'^^ given by the n = N — 2 equation of ( |A.3|) , 



(A.3) 



R 



(2) 
N 



-iq 



1 — a 



7,-N 



[I + - e*«(v^)-Vjv(2Si?!i^ - e 



N-2 



m=l 



(A.4) 



n-l 



The Lyapunov exponent is defined by (with (fN-n = Y\_ RN-mQi ^ 



"Pi. 



m=0 



A = lim — In \(pN- 

{n<N)^oo nC 



(n<Af)^oo nC > ^ — ^ 



(A.5) 



where \n\RN-p\ is expanded, as in ([ToD , to second order in the site energies and 
the term proportional to ln(|Qi|) vanishes for n ^ oo. The disorder averages 



{R^p), {Rn'''), (Rn'-p) and {R^j^>) involved in the expansion of (ln|i?7v-p|) are 
computed from the exact solutions (|A.1| - |A.4| ) and are found to coincide with the 
site-independent values (20. a, b), as expected. From these results we then obtain 



?(1)2\ 



?{2) 



?(2) 



A 



^0 « 



4 (l-a)2 \ 1-5-^ 



1 + a 



+ C.C. 



(A.6) 



which coincides with minus the Lyapunov exponent (|2lD along the opposite direction 
on the ring, as expected. 
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